Module 22: Simple Harmonic Oscillation and Torque 

22.1 Introduction 

We have already used Newton's Second Law or Conservation of Energy to analyze 
systems like the bloc-spring system that oscillate. We shall now use torque and the 
rotational equation of motion to study oscillating systems like pendulums or torsional 
springs. 

22.2 Simple Pendulum 

A pendulum consists of an object hanging from the end of a string or rigid rod pivoted 
about the point S . The object is pulled to one side and allowed to oscillate. If the object 
has negligible size and the string or rod is massless, then the pendulum is called a simple 
pendulum. The force diagram for the simple pendulum is shown in Figure 22.1. 




Figure 22.1 A simple pendulum. 

The string or rod exerts no torque about the pivot point S . The weight of the object has 
radial f - and 6 - components given by 

mg = mg(cosdr -sin09) (22.2.1) 

and the torque about the pivot point S is given by 

x s = f s m x mg = I r x mg(cosd f - sin06) = -Imgsind k (22.2.2) 

and so the component of the torque in the z -direction (into the page in Figure 22.1 for 6 
positive, out of the page for 8 negative) is 

(t s \ = -mglsmd . (22.2.3) 
The moment of inertia of a point mass about the pivot point S is 
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I s = ml 2 . 



(22.2.4) 



From the rotational dynamical equation is 

r ^ J J d2 ° 
( T .<X = I, a = /, 



S'z S S 1.2 

at (22.2.5) 

, • a i2 d 2 e 
-me/ suit? = ml — -. 

dt 2 

Thus we have the equation of motion for the simple pendulum, 

l|. = _£ s in0. (22.2.6) 
dt I 

When the angle of oscillation is small, then we can use the small angle approximation 

sin0 = 0; (22.2.7) 
the rotational dynamical equation for the pendulum becomes 

(22.2.8) 

dt 2 I 

This equation is similar to the object-spring simple harmonic oscillator differential 
equation from (add reference), 

d jc k s~.~. — . 

— T = x, (22.2.9) 

dt m 

which describes the oscillation of a mass about the equilibrium point of a spring. Recall 
that in (add reference), the angular frequency of oscillation was given by 



% ri n g -J-- (22.2.10) 
By comparison, the frequency of oscillation for the pendulum is approximately 

% e „ dulum -^f , (22.2.11) 

with period 
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lit I 

T = = 2jt -. (22.2.12) 



^pendulum 



A procedure for determining the period for larger angles is given in Appendix 22 .A. 
22.3 Physical Pendulum 

A physical pendulum consists of a rigid body that undergoes fixed axis rotation about a 
fixed point S (Figure 22.2). The gravitational force acts at the center of mass of the 
physical pendulum. Suppose the center of mass is a distance / cm from the pivot points* . 




Figure 22.2 Physical pendulum. 

The analysis is nearly identical to the simple pendulum. The torque about the pivot point 
is given by 

T„ = ?„ x fflg = / f x me(cos0 f - smd 6) = -/ mesin^k. (22.3.1) 

5 S, cm ° cm ° v ' cm ° v ' 

Following the same steps that led from Equation (22.2.2) to Equation (22.2.6), the 
rotational dynamical equation for the physical pendulum is 

( \ j j d2 ° 

ai (22.3.2) 

Thus we have the equation of motion for the physical pendulum, 
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(22.3.3) 



As with the simple pendulum, for small angles sind « 6 and Equation (22.3.3) reduces to 
the simple harmonic oscillator equation with angular frequency 



It is sometimes convenient to express the moment of inertia about the pivot point in terms 
of l cm and I cm using the parallel axis theorem (add link), with d S cm = l cm , 

7 5 = 7 cm + ml L , with the result 



Thus, if the object is "small" in the sense that 7 cm <sc ml 2 cm , the expressions for the 

physical pendulum reduce to those for the simple pendulum. Note that this is not the case 
shown in Figure 22. 1 1 . 

22.3.1 Example: Oscillating rod 

A physical pendulum consists of a uniform rod of length d and mass m pivoted at one 
end. The pendulum is initially displaced to one side by a small angle 6 and released 

from rest. You can then approximate sin0 = 6 (with measured in radians). Find the 
period of the pendulum. 




(22.3.4) 



and period 




(22.3.5) 




(22.3.6) 




We shall find the period of the pendulum using two different methods. 



1 . Applying the torque equation about the pivot point. 
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2. Applying the energy equation. 
Applying the torque equation about the pivot point. 

With our choice of rotational coordinate system the angular acceleration is given by 

d = ^k. (22.3.7) 
dt 

The force diagram on the pendulum is shown below. In particular, there is an unknown 
pivot force, the gravitational force acting at the center of mass of the rod. 




f 

The torque about the pivot point is given by 

^ = r ,, em x™g • (22-3.8) 

The rod is uniform, therefore the center of mass is a distance d 1 2 from the pivot point. 
The gravitational force acts at the center of mass, so the torque about the pivot point is 
given by 

% p = (d/ 2)f x mg(-sind 8 + cosf ) = -(d 1 2)mgsind k . (22.3.9) 
The rotational dynamical equation (torque equation) is 

x p = I p a. (22.3.10) 

Therefore 

d 2 d 

-(d/2)mgsmdk = I p —k. (22.3.11) 

dt 

When the angle of oscillation is small, then we can use the small angle approximation 
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Then the torque equation becomes 



sin0 = 6 . 



(22.3.12) 
(22.3.13) 



which is the simple harmonic oscillator equation. The angular frequency of oscillation 
for the pendulum is approximately 



-."f^- (22.3.14) 

The moment of inertia of a rod about the end point P is I p = (1 / 3)md 2 therefore the 
angular frequency is 



(^1(3/2)* 
V(l/3)mJ 2 V <i 



with period 



T- — m2jtJ--. (22.3.16) 



Applying the energy equation. 



Take the zero point of gravitational potential energy to be the point where the center of 
mass of the pendulum is at its lowest point, that is, = 0. 




When the pendulum is at an angle 6 the potential energy is 
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U = mg^(l-cos0). (22.3.17) 



The kinetic energy of rotation about the pivot point is 



2 

The mechanical energy is then 



K rot =-I p co 2 . (22.3.18) 



£ = C/ + ^ rot =mg|(l-cos0) + i/y, (22.3.19) 

with I p = (l/3)mJ 2 . There are no non-conservative forces acting, so the mechanical 
energy is constant, and therefore its time derivative is zero 

n dE d . , dd T dco ,„„ „ „ n . 

= — = me— sine — + Ioj . (22.3.20) 

dt 2 dt p dt 

Recall that co = dd I dt and dco I dt = d 2 6 1 dt 2 so Eq. (22.3.20) becomes 



d ■ n r d 2 d 
= mg — smda> + Icq — -. (22.3.21) 

2 " dt 2 

Therefore two solutions, <x> = , in which the same remains at the bottom of the swing 
and 

O = mg-sin0 + /„^-f (22.3.22) 

2 " dt 2 

Using the small angle approximation, we have the simple harmonic oscillator equation 
(Eq. (22.3.13)) 

^ mg{ dii ) g _ _ (22323) 

dt I p 

22.3.2 Example: Torsional Oscillator Solution 

A disk with moment of inertia about the center of mass I cm rotates in a horizontal plane. 
It is suspended by a thin, massless rod. If the disk is rotated away from its equilibrium 
position by an angle , the rod exerts a restoring torque given by r cm = -yd .At t = , 
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the disk is released from rest at an angular displacement of 6 . Find the subsequent time 
dependence of the angular displacement d(t) . 
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Solution: Choose a coordinate system such that k is pointing upwards, then the angular 
acceleration is given by 



a = — -k, 

dt 2 



(22.3.24) 



The torque about the center of mass is given in the statement of the problem as a 
restoring torque 



x = -yd k . 

cm > 



(22.3.25) 



Therefore the k -component of the torque equation x = / d is 

r 11 cm cm 



-yd = I 



d 2 e 



cm dt 1 



(22.3.26) 



This is a simple harmonic oscillator equation with solution 



6(t) = Acos(co 1) + i?sin(a; 1) 



(22.3.27) 



where the angular frequency of oscillation is given by 



to, 



(22.3.28) 



The z -component of the angular velocity is given by 



dd 
dt 



(t) = -co Asm(co Q t) + co Q Bcos(co 1) . 



(22.3.29) 
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The initial conditions at t = , are that 6(t = 0) = A = 6 , and (dQ I dt)(t = 0) = co Q B = . 
Therefore 



Appendix 22.A: Higher-Order Corrections to the Period for Larger 
Amplitudes of a Simple Pendulum 

In Section 22.2 we found that using the small angle approximation the period for a simple 
pendulum is (22.2.12) 



How good is this approximation? If the pendulum is pulled out to an initial angle 6 that 

is not small (such that our first approximation sin0s0 no longer holds) then our 
expression for the period is no longer valid. We then would like to calculate the first- 
order (or higher-order) correction to the period of the pendulum. 

Let's first consider the mechanical energy, a conserved quantity in this system. Choose 
an initial state when the pendulum is released from rest at an angle 6 ; this need not be at 

time t = 0, and in fact later in this derivation we'll see that it's inconvenient to choose 
this position to be at t = . Choose for the final state the position and velocity of the bob 
at an arbitrary time t . Choose the zero point for the potential energy to be at the bottom 
of the bob's swing (Figure 22.A.1). 




(22.3.30) 






Figure 22.A.1 Energy states for a simple pendulum. 



The initial mechanical energy is 



E Q = K + U = mgl(l-cosd ) . 



(22.A.1) 



The tangential velocity of the bob at an arbitrary time t is given by 
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7 dd , 
v tan =/— , (22.A.2) 



and the kinetic energy at the final state is 



2 



1 2 1 ( 1 dO\ 
K f =-mvl=-m\l—\ . (22.A.3) 



dt 



The final mechanical energy is then 



1 / ,J0\ 2 



£ 7 =K f +U f = 2 m [ l ~^ j + 1(1- cosd) . (22.A.4) 

Since the tension in the string is always perpendicular to the displacement of the bob, the 
tension does no work and mechanical energy is conserved, E f = E . Thus 

1 ( dd\ 2 

— m\l — + mgl{\ - cos0) = mgl{\ - cos0 o ) 

2 \ dt ) 



(22.A.5) 



( l ~dtj =2 y( cos0_cos0 o)- 



We can solve Equation (22.A.5) for the angular velocity as a function of 6 , 

^^/cos0-cos0 o . (22.A.6) 



dd_^ \2g 
dt 



Note that we have taken the positive square root, implying that dd /dt>0 . This clearly 
cannot always be the case, and we should change the sign of the square root every time 
the pendulum's direction of motion changes. For our purposes, this is not an issue. If we 
wished to find an explicit form for either 6(t) or t(d) , we would have to consider the 

signs in Equation (22.A.6) more carefully. 

Before proceeding, it's worth considering the difference between Equation (22.A.6) and 
the equation for the simple pendulum in the simple harmonic oscillator limit, 



In both Equations (22.A.6) and (22.A.7) the last term in the square root is proportional to 
the difference between the initial potential energy and the final potential energy. The 
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final potential energy for the two cases is plotted in Figures 22.A.2 below for -it <6 <n 
on the left and -Jt I2<6 <Jt 12 on the right (the vertical scale is in units of mgl ). 
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Figures 22.A.2 Potential energies as a function of displacement angle. 

It would seem to be to our advantage to express the potential energy for an arbitrary 
displacement of the pendulum as the difference between two squares. This is done by 
recalling the trigonometric identity 



l-cos0 = 2sin 2 (0/2) 



(22.A.8) 



with the result that Equation (22.A.6) may be re-expressed as 



^ = J^2(sin 2 (0 o / 2) - sin 2 (0 / 2)) 



(22.A.9) 



(Note that using Equation (22.A.8) is "undoing" one step of (22 A.5).) 
Equation (22.A.9) is separable, 



dd 



= 2A-dt 



yjsin 2 (0 o /2)-sin 2 (6/ 2) V / 



(22. A. 10) 



What comes next may seem like it's been pulled out of a hat. The motivation is 
analogous to for the simple harmonic oscillator. That is, rewrite Equation (22. A. 10) as 



dd 



sin(0 o /2) 1-- 



sin 



sin 



2 (0/2) 
(0„/2) 



= 2 A %dt. 



(22.A.11) 



12/28/2010 11 



The ratio sin(0 12)1 sin(0 o / 2) in the square root in the denominator will oscillate (but 

not with simple harmonic motion) between -1 and 1 , and so we will make the 
identification 

sin^ Sm(g/2) . (22.A.12) 
sin(0 o /2) 

If (p were proportional to the time t , the motion would be simple harmonic motion. We 

don't expect this to be the case, but we might make some simplifications in doing the 
needed integration, and allow comparison to simple harmonic motion. 



Let k = sin(0 o / 2) , so that 



sin — = ksind) 
2 



1/2 (22.A.13) 



^ ( 1 ■ 2@\ /i ,2-2 j\1/2 

cos— = |^1 -sin — j ={\-k sin <p)' . 

Please be sure to note that here k = sm(d 12) is a dimensionless parameter of the system, 
not a spring constant. 

The integral in (22. A. 1 1) can then be rewritten as 

f . d ° =2[\^dt. (22.A.14) 
J k-sj\-sm 2 (j) J V/ 

From differentiating the first expression in Equation (22.A.13), we have that 

1 

—cos— dd =kcosd) dd) 

2 2 

cos</> „, J\-sm 2 (h 

dd = 2k —d(p = 2k / v d<p (22.A.15) 

cos(0/2) ^/l- sin 2 (0/2) 

= 2k J 1 -™ 2 * dtp. 
^\-k 2 sm 2 (p 



Substituting the last equation in (22.A.15) into the left-hand side of the integral in 
(22. A. 14) yields 
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2k 



.ux Jl-sin 2 ^ j p d(p 

I ; 1 d4> = 21 . 

J k yjl - sin 2 $ yjl - k 2 sin 2 <j> J yjl-k 2 sin 2 <j> 



(22.A.16) 



Thus the integral in Equation (22.A.14) becomes 



yj\-k 2 sm 2 (f) 



g 



-w 



dt 



(22.A.17) 



This integral is one of a class of integrals known as elliptic integrals. We will encounter a 
similar integral when we solve the Kepler Problem, where the orbits of an object under 
the influence of an inverse square gravitational force are determined. 

We find a power series solution to this integral by expanding the function 



(1 - k 2 sin 2 0)- 1/2 = 1 + - k 2 sin 2 <j> + - k 4 sin 4 </> + 
2 8 



(22A.18) 



The integral in Equation (22A.17) then becomes 



jj 1 + hi 2 sin 2 <j> + ^k 4 sin 4 + • • • j dcp = dt . (22 A. 1 9) 



Now let's integrate over one period. Set t = when d = , the lowest point of the swing, 
so that sin0 = and (p = . One period T has elapsed the second time the bob returns to 
the lowest point, or when (j) = 2it . Putting in the limits of the (j) -integral, we can 
integrate term by term, noting that 



r?- n 1 



1 . 2 1 



J q - k 2 sin 2 (p d<j) = J q -£ 2 -(l-cos(20))rf0 



sin(2^ 
2 l( 2 ) 

1 , 2 1 • 2 ^0 

= —Jtk = — jrsin — . 

2 2 2 



2 ST 



(22.A.20) 



Thus, from Equation (22.A.19) we have that 
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1 



2tc + —jt sin 2 — + • • • = 



(22.A.21) 



We can now solve for the period, 




(22.A.22) 



If the initial angle is small compared to 1 (measured in radians), d « 1 , then 
sin 2 (0 o 12) = d 2 1 A and the period is approximately 



T*2a.l-(\ + ^\ = T 



where 



/ 1 7 \ 
1 + — 6/ 2 



16 



(22.A.23) 



T =2a. 



(22.A.24) 



is the period of the simple pendulum with the standard small angle approximation. 
The first order correction to the period of the pendulum is then 



(22.A.25) 



Figure 22 .A. 3 below shows the three functions given in Equation (22.A.24) (the 
horizontal, or red plot if seen in color), Equation (22.A.23) (the middle, parabolic or 
green plot) and the numerically-integrated function obtained by integrating the expression 
in Equation (22.A.17) (the upper, or blue plot) between (p = and (j) = 2it as a function 

of k = sin(0 o / 2) . The plots demonstrate that Equation (22.A.24) is a valid 

approximation for small values of d , and that Equation (22.A.23) is a very good 

approximation for all but the largest amplitudes of oscillation. The vertical axis is in 

units of yfT/ g . Note the displacement of the horizontal axis. 
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Figure 22.A.3 Pendulum Period Approximations as Functions of Amplitude. 
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